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The conjugate pairing of Lyapunov exponents for a field-driven system with smooth inter- 
particle interaction at constant total kinetic energy was first proved by Dettmann and Mor- 
riss [Phys. Rev. E 53, R5545 (1996)] using simple methods of geometry. Their proof was 
04 ■ extended to systems interacting via hard-core inter-particle potentials by Wojtkowski and 

(«_^ ' Liverani [Comm. Math. Phys. 194, 47 (1998)], using more sophisticated methods. Another, 

>--^ , and somewhat more direct version of the proof for hard-sphere systems has been provided by 

^^ ' Ruelle [J. Stat. Phys. 95, 393 (1999)]. However, these approaches for hard-sphere systems 

P^ I are somewhat difficult to follow. In this paper, a proof of the pairing of Lyapunov exponents 

^ ■ for hard-sphere systems at constant kinetic energy is presented, based on a very simple ex- 

'~~i ' plicit geometric construction, similar to that of Ruelle. Generalizations of this construction 

^^3 , to higher dimensions and arbitrary shapes of scatterers or particles are trivial. This construc- 

CN| ' tion also works for hard-sphere systems in an external field with a Nose-Hoover thermostat. 

However, there are situations of physical interest, where these proofs of conjugate pairing 
^m^ , rule for systems interacting via hard-core inter-particle potentials break down. 
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^ '. I. INTRODUCTION 

^^ ■ Thermostattcd, field-driven systems have been popular models for non-equilibrium molecular dynamics 

^ I (NEMD) simulation studies of transport processes in fluids. NEMD studies consider systems with a large 

OO ■ number of particles interacting with each other, driven by an external field ||l|,pi. In these studies, the 

thermostat continuously removes the energy generated inside the system due to the work done on it by the 
external field, by means of a dynamical friction term in the equations of motion. One finds that for these 
f~^ I systems, a non-equilibrium steady state (NESS), homogeneous in space, is reached after a sufficiently long 

^ ■ time |H-p[. Here we consider a particular kind of thermostat, where the friction is linearly coupled to the 

^ I laboratory momenta and keeps the total laboratory kinetic energy of the particles constant (i.e, isokinetic 

^ ■ Gaussian thermostat coupled to the laboratory momenta). The dynamical description of these systems, would 

'~| I be Hamiltonian in the absence of this thermostat, and the corresponding Hamiltonian system, obtained by 

^ ■ dropping the dynamical friction term from the equations of motion, will be referred to as the background 

^ I Hamiltonian system. Due to the presence of the dynamical friction terms in the equations of motion, such 

a system is no longer Hamiltonian with a conservation of phase space volumes, but instead is phase-space 
contracting |0,p[. The sum of all the Lyapunov exponents, which measures the rate of long-time exponential 
i^ ■ growth of phase-space volume, for such a system is thus negative l^. Of course, the background Hamiltonian 

system, if chaotic, is such that the Lyapunov exponents sum up to zero. Furthermore, due to its symplectic 
form, the Lyapunov exponents of the background system come in pairs such that the sum of each such pair 
is also zero |10|T^. The phenomenon of such pairing of Lyapunov exponents, where the sum of each pair 



o 



X 



of non-zero exponents takes a constant value independent of the particular pair, is known as the conjugate 
pairing rule (CPR). Dettmann and Morriss p2| have studied the isokinetic Gaussian thermostattcd field- 
driven systems that are under consideration here, assuming that the particles of the system interact with 
smooth pair potential energies, and that the forces on the particles due to the external field depends only on 
their positions. They have proved that under these conditions, in a restricted subspace of the phase space of 
all the particles, characterized by all the non-zero Lyapunov exponents, such a system is /i-symplecticFl As a 



^For a definition of /.t-symplecticity condition, see Eqs. (w) and (Mj of this paper. The usual symplectic condition is 
a particular case of Eq. (H) with fi = 1. 



consequence, the CPR is exactly satisfied in that subspace, and it is independent of the number of particles 
in the system (corresponding simulation results can be found in [ ^3^141 ]) — the sum of each pair comes out to 
be the same negative constant. One important consequence of this result is that the macroscopic transport 
coefficients of these systems, in the linear order, can be obtained from this constant value of the sum (see, for 
example, |1^]). The restricted subspace, characterized by all the non-zero Lyapunov exponents, is identified 
by observing that with an isokinetic Gaussian thermostat, trajectories of the system always lie on a constant 
total kinetic energy hypersurface in the phase space of all the particles. This constraint generates a zero 
Lyapunov exponent. Also, two points in the phase space do not separate exponentially in time if one point 
is chosen in the direction of flow of the other. This generates another zero Lyapunov exponent. 

As hard-core inter-particle interaction can be dealt as a limiting case of a very short range smooth potential, 
one would expect that in this restricted subspace, the system is still /it-symplectic, and as a result, the CPR 
will continue to hold for a field-driven isokinetic Gaussian thermostatted system with hard-core inter-particle 
interactions. The corresponding pairing of Lyapunov exponents for hard-sphere systems has been proved 
by direct means using the differential geometric structure of the phase space, by Wojtkowski and Liverani 
|l6| . Another version of the proof for hard-sphere systems has been obtained by Ruelle jl^l . The above two 
approaches for hard-sphere systems are somewhat difficult to follow. The purpose of this paper therefore, 
is to present a proof of the pairing of Lyapunov exponents (also by direct means) for hard-sphere systems 
at constant kinetic energy, based on a very simple explicit geometric construction, similar to that used by 
Ruelle. Two kinds of field-driven isokinetic Gaussian thermostatted systems with hard-core interaction are 
considered here: in section II, the proof of the CPR is carried out for the three-dimensional Lorentz gas, 
where mutually non-interacting point particles suffer specular collisions with fixed spherical scatterers. In 
section III, the proof is then carried out for a gas of hard spheres. The explicit method of this paper allows 
one to identify the dependence of these approaches (described in Refs. |l^,|l3l and here) on the geometrical 
shapes of the scatterers (or the particles, as the case may be) and on the nature of the externally applied field. 
Based on it, generalizations to higher dimensions, to arbitrary geometry of the scatterers or the particles, 
and to the case where the masses of the particles are arbitrary, become trivial. Finally, in Sec. IV, we 
argue that this construction can be used to prove the CPR for hard-sphere systems in an external field with 
a Nose- Hoover thermostat. We also identify situations of physical interest, where these approaches break 
down. 



II. PROOF OF THE CPR FOR THREE-DIMENSIONAL LORENTZ GAS 

The Lorentz gas model consists of a set of scatterers fixed in space together with mutually non-interacting 
moving particles that suffer elastic, specular collisions with the scatterers. Here we consider the version 
of the model in three dimensions where the scatterers are hard spheres, and are placed in space without 
overlapping. Each of the moving particles is a point particle with unit mass (m = 1), and is subjected to 
an external force that depends only on its position, as well as a Gaussian thermostat which is designed to 
keep its kinetic energy at a constant value. During a flight, the equation of motion of a particle is 



r, p 



p,F— ap 



(1) 



where F = ~ S/(j) is the force on the particle due to the external field and a is the coefficient of dynamical 

friction representing the isokinetic Gaussian thermostat. The value of a is obtained from the fact that the 

P^ . 
kinetic energy of the particle — is constant during a flight, i.e, 

« - ^ (2) 

and without any loss of generality, the time hereafter is rescaled such that p^ = 1. At a collision with 
a scatterer, the post-collisional position and momentum of a particle, F^ — [r^, p_|_], are related to its 
pre-collisional position and momentum F_ = [r_, p_], by 

F+ = QF_ == [f_, p_ - 2(p_ -fi)?!] , (3) 

where n is the unit vector from the center of the scatterer to the point of collision (see Fig. 1). 



It is sufBcient to consider the motion of only one of the moving particles in its six-dimensional phase space 
to investigate the chaotic properties of this system, since they do not interact with each other. To obtain the 
Lyapunov exponents, we consider a particle at the phase space location Fq = [fo, po] at time i = 0. In time 
t, it suffers s sequential collisions at time instants ti,t2,--- ,ts with the scatterers. In between collisions, 
the particle undergoes flights, acted upon by the external field. We refer to the phase space trajectory of 
it as the "reference trajectory". We also consider another particle at Fq + STq = [tq + Sro, po + Spg] at 
time i = 0, such that Fq and Fq + STq are infinitesimally apart from each other. This particle suffers the 
same sequence of collisions and likewise its movement in the phase space generates the "adjacent trajectory" . 
The two trajectories remain infinitesimally apart from each other at all times. Let H(tj — tj~i) denote the 
time evolution operator for ST{t) due to a flight between time instants ij_i and tj and let M^ denote the 
evolution operator for 6r{t) at the i-th collision. For < ii < ^2 < ■ ■ • < ^s < i, we therefore have 

ST{t) = H(t - ts) Ms Hits - ts-i) Ml H(ii) <5ro = L(i) ^Fq . (4) 

Let us now also define another six-dimensional matrix T(t), such that during a flight 

5T{t) = T(i) 6T{t) . (5) 

The matrix L(t) can then be obtained from the solution of the differential equation 

H(t) = T(t)H(i), (6) 

and Eq. (Q), with the boundary condition that H(0) = L(0) = I. The Lyapunov exponents, measuring the 
rate of exponential separation between the reference point and adjacent point in this six-dimensional phase 
space for long times, are then defined as the logarithms of the eigenvalues of the matrix A, where 

A= lim{[L(0FL(i)}V2*. (7) 

t — >oo 

Clearly, there can be at most six Lyapunov exponents of this system, two of which are zero due to the 
reasons explained in the introduction. For our purpose, we need to select out a four-dimensional subspace of 
this six-dimensional phase space where all the Lyapunov exponents may be non-zero. In this four-dimensional 
subspace, the proof of the CPR would follow from the /i-symplecticity property of L(i) |12,nq,0]. Thus, all 
we need to prove to establish the CPR in this four-dimensional subspace is that there exists a /i(i), such 
that with the four-dimensional subspace representation of L(t), 

A*(t)[L(t)]Tj[L(t)] = J, (8) 

which is the /i-symplecticity condition. Here, J is a 4 x 4 matrix with each entry being a 2 x 2 matrix : 



Our goal is to prove Eq. (ph by means of the geometric construction, developed in Rcf. |12|. 

The identification of this four-dimensional subspace is facilitated by decomposing the six-dimensional 
position and momentum space of the particle at every point on the reference trajectory, into two separate 
three-dimensional subspaces, one for the position-space and the other for the momentum-space. Next, in 
each of these two three-dimensional subspaces, a unit vector eo = p is chosen and two other unit vectors 
ei and 62 are also chosen to form a complete set of orthonormal basis. One zero Lyapunov exponent, which 
occurs due to the fact that 5p must be chosen orthogonal to eo in order to respect the constraint that p^ = 1, 
is avoided by measuring 5p by its components along the local directions of ei and 62, i.e, 

2 
Sp ^ ^ Spi ii . (10) 

1 = 1 

The other zero Lyapunov exponent, which occurs due to the fact that the adjacent point does not exponen- 
tially separate from the reference point if Stq is chosen along eo(i = 0), is avoided by measuring dr also by 
its components along the local directions of ei and 62, i.e. 



Sr — y^ Sr.i I 



(11) 



Albeit Co, being the momentum of the particle, is uniquely defined at each point on the reference trajectory, 
ei and 62 can be chosen arbitrarily at every point of the reference trajectory, maintaining the orthonormality 
condition. This ambiguity in the local orientations of ei and 62 can be removed by an initial choice of ei (t — 0) 
and 62 (t = 0) at {rQ,po) and subsequently connecting the orthonormal set of basis vectors at different points 
on the reference trajectory by means of a "parallel transport". The construction and parallel transport of 
these basis vectors are the key components of the proof of the /z-symplecticity in this procedure. 

To this end, we first concentrate on the evolution of 6r{t) during the process of a flight of the particle, 
which is a special case of the systems that Dettmann and Morriss [|l2| have considered, where (j) is the 
potential due to only the external field. Following their construction, we therefore use 



e* = - {F ■ e.i)eo, i = 1,2 

while the parallel transport of cq is obtained from the equations of motion, Eq. (m), i.e, 

2 



eo = ^(F ■ ^i) ^i 



eo 



(12) 



(13) 



From Eq. Im), having obtained 



eo , 



it is then easy to show that during a flight 



E(^^ 



(14) 



Spt and Sp.^ =^ [- ViV jcj) - {F ■ ei){F ■ Cj)] Srj - a Spi 



(15) 



Consequently, the matrix T can be expressed as 



I 



(16) 



where each of the elements in T above is a 2 x 2 matrix and R/ is a symmetric matrix. The matrix T in 
Eq. ( p^ has the property that 



T^J + JT = aJ, 

from which it can be easily shown that for a flight between tj _ 1 and tj 



li(t, - i,-i)[H(i, - i,_i)]Tj[H(t, 



h-i) 



(17) 



(18) 



where /^(tj — tj-i) — exp j^' a{t')dt' 



along the reference trajectory. 

Next, we study the evolution of the infinitesimal volume element 5T due to a collision with a scatterer to 
obtain a similar four-dimensional representation of the matrix M, defined in Eq. (0). We notice that even 
though Eqs. (O) and (ITsI) uniquely determine the orientations of the basis vectors over a flight given their 
orientation at the initiation of the flight, the equations of parallel transport connecting the pre-coUisional 
and post-coUisional basis vectors arc still lacking. As the post-collisional momentum of the particle eo-j- can 
be related to its pre-collisional momentum eo_ by 



eo+ 



eo- - 2(eo- • n) n , 



(19) 



over a collision [see Eq. (^)], a parallel transport of the orthonormal basis vectors that serves our purpose, 
can be completed also by using 



-i+ 



2 (ei_ ■ fijfi . 



1, 2 



(20) 



With the use of Eqs. ( |lO| - |ll| ) and ( |l9| - P0| ), wc now show that as the infinitesimal pre-coUisional phase space 
separation can be written as 



5V. 



6r- , (5p_ 



y^ Sn- Ci- , ^ Spi- ii 



. i = 1 i = 1 

the corresponding infinitesimal post-coUisional separation, (5r_|_, can then be expressed as 



ST. 



Sr+ , Sp^ 



2 2 

)ri+ ei+ , ^ (5pi+ ei+ 

4=1 1=1 



E^^ 



(21) 



(22) 



which would also reduce the collision dynamics oiST{t) to four dimensions. We want to obtain the symplectic 
property of the 4x4 matrix M, defined by 



<5r. 



M5r_ 



(23) 



We begin by studying Fig. 2, which is the two-dimensional projection of an exaggerated picture of a 
collision that is taking place in three-dimensions. The plane A_, perpendicular to the reference trajectory 
at its point of collision A intersects the adjacent trajectory at B. While (5r_ — AB, Sp_ too lies on the plane 
A_. Similarly, A+ is the plane that is perpendicular to the reference trajectory at D passing through C, the 
point of collision of the adjacent trajectory. The post-coUisional position-space separation between the two 
trajectories, Sr+ = DC and 5p_^_ also lies on the plane A+. Clearly, there is a time gap St between the two 
collisions at A and C (i.e, the time required for the adjacent point to travel from B to C), given by 



6t 



Sr^ ■ h 



(24) 



Following the procedure outlined in |Q , we find that the infinitesimal phase space separation of the two 
trajectories just before the two collisions at A and C is 



ST* = 5r_ 



r_(5T 



and consequently. 



5T, 



dQ_ 

(9r_ 



5r* 



T+St 



(25) 



(26) 



where r_ (r+) is obtained from the equations of motion of the particle for the flight immediately before and 
after the collision at A [see Eq. (|lj)], i.e. 



eo±, ^{F* ■ei±)e,± 



{F* is the force on the reference point due to the external field at the point of collision at A) and 



(27) 



dQ_ 

9r_ 



2eo-- 



I 




1 


' dfi ^ ^ dn 
df- dr- 


I 


- 2nn 



(28) 



[obtained from Eq. (|^)]. Starting with the expression of 5r+, we obtain, from Eqs. (|lj) and ( p4[]2q ) that 

(5r+ = (5r_ - 2 ((5r_ • n) n , (29) 

which, together with Eqs. (|20|-p2[) yields 

eo+ • 6r+ = and 6rij^ = 5ri- . i — 1, 2 (30) 



In a similar manner, the expression of 5pj^ can also be easily obtained from Eq. (|26|-p8|) , i.e, 

2 2 

5p+ = A- 5r*_ + (I - 2hn) Sp_ + (I - 2hh) ^(F* • e,_) ei_ St - ^l^^* ' ^*+) ^■'+ "^^ ' (^1) 



j = i 



where A = — 2eo- 






and Sr_ = {Sr 



eo- St) is the infinitesimal vector AC lying on 



the surface of the scatterer (see Fig. 2). 

The simplification of the expression on the r.h.s. of Eq. ( pi| ) is carried out in Appendix A. After collecting 
all the terms together from Eqs. (Al), (A2) and (A5), we obtain 



Spi- 



E 



(5y Spi^ 



{F* ■ h)ie^^ ■n)iej^ ■ n) 1 (e^- ■ n)(ej_ ■ fi) + Sj^j (ep- ■ n)^ 

(eo- ■ fi) a (eo_ ■ n) f J- 



(32) 



One can now see from Eqs. (E3|), (29) and (|32|) that the 4x4 matrix M has the following structure 



M 



I 



(33) 



where the 2x2 matrix Re is symmetric. This allows us to conclude that in the four-dimensional subspace 
representation, the transformation of ST over a collision is symplectic, i.e, 

M'^JM = J. (34) 

Finally, using Eqs. ([^, (|lq), (|34|) and ^{t) = exp La{t')dt' , the four-dimensional subspace representation, 

L(i) is seen to be /i-symplectic. The four non-zero Lyapunov exponents of this system can therefore be 
arranged in pairs such that the sum of each pair is exactly equal to — (a)t, the negative of the long time 
average of a along the reference trajectory, i.e., the CPR is exactly satisfied for this system in the restricted 
four-dimensional subspace. 



A. Further Generalizations of this Geometric Construction 



We end this section with the following generalizations: 
(i) The proof can be extended to any arbitrary dimensions, simply by including appropriate number of unit 
vectors e^'s. 

(a) Equation (pi) demonstrates the role of geometry of the scatterers and the effect of the external field 
on the symmetry property of Re, as having a symmetric Re is essential for the symplecticity of M [and 
consequently, the /i-symplecticity of L(i)]. The effect of the geometry of the scatterers is coded in the 
A • Sr_ term of Eq. (^l|). For nonspherical scatterers, this geometric construction works, so long as the 
surfaces of the scatterers are smooth. One has to notice that at the collision point A (see Fig. 1) on an (ti — 1)- 
dimensional surface embedded in an n-dimensional Euclidean manifold, Sr_ is the infinitesimal vector AC 



along the surface of the scatterer, and since Sn appearing in Eq. (A4) is the infinitesimal difference between 
the two unit normal vectors to the surface of the scatterer at A and C, for a smooth surface, one can define 
an n X n matrix B, such that Sn = 'B ■ Sr_. In general, the form of the matrix B depends on the shape of 
the surface of the scatterer at collision point A. As a special case, if the scatterer is spherical, then B can be 
explicitly constructed to be the identity matrix times the inverse radius of curvature of the sphere. In terms 
of B, one then simply needs to obtain the form of A • Sr_, analogous to Eq. (|A5|), given below as 



A-Sr_ ^ 2 Y^ 

i,3 = i 



(ej_ • n)(ej_ • B • cq-) + (e^- • n)(ei_ • B • cq-) 

- (eo- • n)(ei- • B • eo_) (e^- • n)(ej- • h) 

eo- ■ n 



)rj^ ej+ 



(35) 



The term in square bracket in Eq. (^5|) is symmetric in i and j, which contributes the symmetry of the 
matrix Re in Eq. (33). However, the role of external field on the symmetry of Re is a bit more subtle — 
the i^*-dependent term in Eq. ( |32| ) arises due to the fact that the reference and the adjacent trajectories 
do not collide the same scatterer at the same instant, despite the fact that collisions between the particles 
and the scatterers are instantaneous. As Eq. ( pTf ) shows, we have used the fact that both the dynamics of 
r_ and r+ involve the same force F* , which is possible only if the pre- and the post-coUisional values of 
F* (i.e., appropriate F* and FJ) arc the same. One can therefore conclude that this construction can be 
generalized to prove CPR for arbitrary (smooth) shapes of the scatterers in any dimensions, when the forces 
on the particles due to the external field depends only the particles' positions. 

(Hi) The proof does not depend on the specific locations of the scatterers, and hence the CPR holds for any 
arrangement of the scatterers in space. 

III. PROOF OF THE CPR FOR HARD SPHERE GAS IN THREE DIMENSIONS 

In this section, we consider a gas of N identical moving spheres in three-dimensions. Each of the spheres 
has a unit mass and is subjected to an external force that depends only on its position, as well as a Gaussian 
thermostat which keeps the total kinetic energy of the system at a constant value. The spheres interact with 
each other by means of binary elastic collisions. During a flight, where no collision takes place between any 
two of the spheres, the equations of motion of the system are 



Pi' 



Pi- 



F' 



apr 



I...N 



(36) 



where Fi> is the external force on the i'-th sphere and a is the coefficient of dynamical friction representing 
the isokinetic Gaussian thermostat (hereafter primed indices will always indicate sphere numbers) . The value 

of a is set in a way that it keeps the total kinetic energy of the system, > — ^ constant, i.e, 

i' = 1 



AT 
Li' = l 



•Pr 



N 
i' = 1 



(37) 



N 



As before, we rescale the time such that > p^, = 1. At a collision between the I'-th and the j'-th spheres. 



the post-collisional positions and momenta are related to the pre-collisional values by 



rj'+ 



Pt'+ 



m- 



Pr-)-ni'j'}fii'y 



and 



Pf+ = Pr- + {{Pi'- - Pf-) ■ hej'}hi,j, , 



(38) 



where fiiiji is the unit vector along the line joining the center of the i'-th sphere to the j'-th sphere at the 
time of the collision. Since such a collision is instantaneous, during any such binary collision process, the 
positions and momenta of the spheres not participating in the collision do not change. 

We define the 3iV-dimensional vectors R and P, assembled respectively from the three-dimensional position 
and momentum vectors of the N spheres, such that 



R ^ [n, r2,. 



rN 



] and P = [pi , p2 



Pn _ 



(39) 



with P~ = 1. We also form two more 3iV-dimensional vectors F and N.i'ir. F describes the external force 



on the spheres, and Ni'j' is assembled from the unit vector n, 



i'j' ; 



Fo 



P 



N 



and N, 



i-y 



1 

V2 



0, 0,. 



(40) 



such that the i'-th and the j'-th entries of Ni'j' (they are the only non-zero entries) are fiiiy j^ji and 
—hiiji/\/2 respectively [satisfying the normalization condition Ni' j' ■ Ni'j' ~ 1]. Equations (pq-BTf) can now 
be rewritten [in the same form as Eqs. (|l|) and (pi)] as 



R, P 



P, F ~ aP 



with 



FP. 



(41) 



while the collision dynamics can be rewritten as 



r+ = Qr_ = 



R- 



P- 



2{P- 



N,'j')N,,,, 



(42) 



analogous to Eq. (H). As discussed in the introduction, the geometric construction associated with Eqs. ( p9| - 
p2| ) can be obtained as an extension of the corresponding constructions for the three-dimensional Lorentz 
gas. We will see that the construction of the unit normal vector Ni'j' is the key component of the proof of 
/^-symplecticity and the CPR for hard-sphere gases. 

The proof of the CPR proceeds exactly in the same way as it has been described in section II. The time 
evolution of the infinitesimal 6A^-dimensional phase space separation T(t) between the reference and the 
adjacent trajectories can be decomposed by means of 6A^-dimensional H and M matrices as in Eq. (|j). 
The (6iV — 2)-dimensional reduction of H matrices can be subsequently obtained by constructing N 3N- 
dimensional basis vectors eo(t = 0), ei(i = 0), ., ., ., ., e(3jv_i)(i = 0) at Fq, and then parallelly transporting 
them using equations analogous to Eqs. (Ojll) and choosing to measure both 6R and SP in directions 
orthogonal to cq in the same manner as in Eqs. (Oxl-O). As a special case of what Dettmann and Morriss 



considered ||13] , in terms of the (6iV — 2)-dimensional representation, the matrix H(ij — tj-i ) is easily seen to 
be /i-symplectic with fi(tj 



tj-i, 



exp 



/j^ a{t')dt' along the reference trajectory, for a flight between 
tj- 1 and tj. What remains to be shown, therefore, is that the matrix M describing the transformation of 
5T(t) due to a binary collision is also symplectic. To this end, following section 11 of this paper, we use a 
simple extension, comprising of the same form of parallel transport (as in Eqs. (|9|) and ([2^)) of the N basis 
vectors cq, ei, ... e(3^_i') over the binary collision between the i'-th and the j'-th sphere, i.e. 



eo+ 



io- ~ 2(eo^-N,,j,)N,, 



and 



=i+ 



2 (e,_ • N,,y) N,,y . l<i< (3iV-l) (43) 



Evaluation of the expression of M in this (6A^ — 2)-dimensional representation is a bit more involved than 
that presented in Sec. II. Nevertheless, one can still use Fig. 2 to illustrate the collisions at A and C and 
the construction of the planes A_ and A-|_ in a schematic way, keeping in mind that this time the diagram 
describes quantities in 3iV-dimensions (for another version of Fig. 2 in the context of hard-sphere gases, the 
reader may also find Fig. 6 in [H helpful). We will follow the logical steps described in Eqs. (p4-28) to 
obtain the expressions of 5X^ and 5P+. Let us first define below the 6iV dimensional vector SY'± describing 
the infinitesimal post(pre)-collisional phase-space separations of the two trajectories 



6T+^ 



SR±, SP± 



Sri± , . . . , SrN± , Sp^^ , . . . , 6pj^^ 



37V-1 



3Ar-l 



y^ SRi± h± , ^ 5Pi± ei± 



(44) 



We observe that the time lag between the binary collisions, involving the i'-ih and the j'-th sphere on the 
reference and the adjacent trajectories (schematically at A and C respectively in Fig. 2) is 



5t 



{Srji^ — 5rifJ) ■ n 



I J 



{pj'- - Pi'-) • hi 



SR- ■ Nij, 
eo_ • Ni'j' 



(45) 



Following the procedure outlined in [[l8| , we find that the infinitesimal phase space separation of the two 
trajectories at A and C is, as before, 



<5r* 



6T. 



r_ St and STj 



dQ_ 

9r_ 



6r* - r+ (5t , 



(46) 



where r_ (r+) is obtained from the equations of motion of the system for the flight immediately before and 
after the collision at A (from Eqs. (|4l|)), i.e. 



r± 



eo± 



E 



{F* ■ ei±)ei± 



(47) 



Here F* is the force on the reference point, described in Eq. 
coUision at A. Using Eq. (B^), 



due to the external field at the point of 



dQ 

9f7 



2eo-- 



■aj^. 



dR- 



m> 



dN,,,, 



OR- 







where each entry of the matrix on the r.h.s. of Eq. (Esh is a 3N x 3iV matrix. 

Starting with the expression of 6R^, once again we obtain, from Eqs. (Bs) and ([45| 



(48) 



that 

6R+ = SR- - 2 {6R- ■ Ni,j,)Ni,j, , (49) 

which, together with Eq. ( [4^ ) yields 

eo+ • SR+ = and SRi+ = 5R,^. l<i< (ZN - 1) (50) 

In a similar manner, the expression of 5P+ can also be easily obtained from Eqs. (^-48), i.e, 
5P+ = A-6R*_ + (I - 2N,>fNi^j')6P- 

+ (I - 2NefN,,j,) Y, {F* ■ e,_) e,_ <5t - ^ {F* ■ e,+ ) e,+ St , (51) 



where A 



2eo_ 



dR- 



■N,, 



N,, 






and 5R_ = {SR- + io-Sr). The simplification of the 



expression of SP+ is carried out in Appendix B. Finally, using Eqs. (B1-B2) and (B6-B7), we obtain that 
in terms of the {6N — 2)-dimcnsional representation, the matrix M can be written in the form 



SP^+ 



3JV-1 

>; 


S^J6P,- + iw,, - 


3JV-1 


i = i 


»,i = i 



{F*-N,,,,)ie,.-Ni,j,){ej.-Ni,j,) 
(eo^ • n) 



6R,- 



(52) 



where Wij = Wji has been evaluated in Eq. (piq) in Appendix B. Equation (p2D, together with Eqs 
implies that M has the same form as in Eq. (|33|), where the {3N — 1) x (37V — 1) matrix Re is symmetric. 
This again allows us to conclude that in the (6iV — 2)-dimensional subspace representation, M is symplectic, 

and L(t) is /i-symplectic with /i(i) — exp L a{t')dt' , i.e., the CPR is exactly satisfied in the restricted 

(6A^ — 2)-dimcnsional subspace for this system. 



A. Further Generalizations of this Geometric Construction 



We end this section with the following observations: 
(i) the condition that each sphere is of unit mass is not essential — if the i'-th particle has a mass rrii' , one 

can define the new quantities f^, = m^,rii, p'-, = m^, '^ pi> and F,^, — m^, ^ Fi' and then begin with the Eqs. 
( Pq - p^ ) replacing fi/'s, p^/'s and F^/'s by the corresponding primed variables. One also needs to use 



N,,, 



0, 0, 



TTii' + rriji 



fc^'^' , . . . , , 



mi' + rriji 







(53) 



such that the z'-th and the j'-ih. entries of N^,-, (they are the only non-zero entries) are 



y 



mil + rrij 



H']' 



fiiiji respectively (satisfying the normalization condition N^,;, ■ N^,,, — 1) and define a 



and — 

y m,ii + TTT-j' 

corresponding (SN x 6N matrix U' (see Eqs. (612) in Appendix B and its preceding paragraph) such that 



rrii' 



I, U' 



'1-3 ^ ^ J « 



(54) 



(ii) the proof can be trivially extended to any arbitrary dimensions by including appropriate number of unit 
vectors e^'s and 

(Hi) Following point (ii) of Sec. II. A, it is easy to generalize this construction, and hence the proof, to 
arbitrary (smooth) shapes of the particles. One simply needs to use an analogous matrix B such that 



6R_. Just from Eq. ([B3D, one then obtains the same equation as Eq. (^ 



A • SR_ = 2 ^ 



(ei_ • Ni,y){ej^ ■ B • eo-) + (ej_ • Nej'){et- ■ B • eo-) 



(eo- • iVi'j')(ei- • B • ej_) -7- (e^^ • Ni,y){ej- ■ Ni,y) 



eo- • iVi'j' 



5Rj- ej+ . (55) 



Notice already the striking similarity of the term inside the square brackets in Eqs. ( pS]) and (55) to the 
form of Wij in Eq. ( B18 ). When the nN x nN matrix B is decomposed into N x N blocks of n x n matrices, 



the only non-zero entries are Bj/^', Hi'j', B^/i' and Bj/^/, corresponding to the two colliding particles i' 
and j' . As a special case, if these particles are spheres, then the matrix B can be explicitly constructed 
to be U/[\/2(ai' + (^j')], where a^' and gji are the radii of the colliding particles [pl| (in Appendix B, we 
have used a^' = Gji = a). The factor of ^/2 follows from Eqs. (^ ) and (B4). Once again, the form of 
Eqs. ( pq ) contributes the symmetry of Re. In a similar manner as explained in point (Hi) of Sec. II. A, the 
symplecticity property of M (and hence the /x-symplecticity of L(t) crucially hinges upon the symmetry of 
Re, for which it is necessary to have the expression of F* to be invariant under a collision. 



IV. DISCUSSION 

While much of the content of this paper is technical, the underlying principle behind our discussion of the 
CPR is similar to those in the existing literature |lq,0]. First, the dimension of the phase space is identified 
where all the Lyapunov exponents are non-zero. The dimension of this reduced phase space, characterized 
by all the non-zero Lyapunov exponents, depends on the number of macroscopic conserved quantities in the 
problem (such as total momentum, total angular momentum, total energy etc.) that are consistent with the 
dynamics. By means of confining the dynamics to a hypersurface of dimension 1 less than that of the full 
phase space, each such conserved quantity reduces the dimension of the phase space characterized by all the 
non-zero Lyapunov exponents by 1. In addition, the fact that two points in the phase space do not separate 
exponentially in time if one point is chosen in the direction of flow of the other, reduces the dimension of the 
phase space characterized by the non-zero Lyapunov exponents also by 1. The dynamics for systems with 
hard-core inter-particle interactions are then decomposed into flights and collisions. The dynamics of a flight 
segment in the reduced phase space is shown to be /i-symplcctic as a special case of the explicit geometric 
construction in Ref. ||l2| . By means of another simple explicit geometric construction, it is also possible to 
incorporate the collisions in this formalism, which subsequently leads one to evaluate the matrix M and 
demonstrate that it is symplectic. Finally, the matrix L(i), assembled from the sequential products of the 
H and M matrices, is easily seen to be /i-symplectic using the fact that the product of a symplectic and a 
/i-symplectic transformation is /i-symplectic. We must note here that the proof of CPR by this construction 
is sufficient and not necessary; and as a result, if a system does not satisfy this proof, it is not guaranteed 
that CPR would not be satisfied for it. Moreover, this explicit geometric construction, used in this paper for 
the parallel transport of basis vectors, is by no means unique. One could conceivably use a different set of 
equations for the parallel transport mechanisms for both flight and collision parts of the dynamics to prove 
the CPR for a dynamical system. 

Questions naturally arise about the applicability of this construction for other hard-particle systems. It 
turns out that using the same explicit geometric construction that is presented here, CPR can also be proved 
for a gas of hard particles (of finite sizes) in an external field (such that the forces on the particles due to 
this field depend only on their positions of the particles) under a Nose-Hoover thermostat. One simply has 
to separate the time evolution of ST in terms of the H and the M matrices in an appropriately reduced 
phase space. The /x-symplecticity of the H matrices can be trivially obtained as a special case of Ref. [EOJ, 
where the CPR has been proved for a system of particles that interact with each other by means of smooth 
interparticle potentials. The symplecticity property of the M matrices can also be easily seen to be valid, if 
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one combines the observation of point (ii) of Sec. III. A with Eq. (15) of Ref. |20|. We have seen in point 
(ii) of Sec. III. A that for M to be symplectic, it is necessary that F* be invariant under a coUision, and it 
is definitely the case with Eq. (15) of Ref. P0| - 

Clearly, this explicit geometric construction allows one to see the role of the external field in between 
collisions (that it is necessary for the force on the particles due to the external field be dependent only on 
their positions), the role of F* and the geometry of the colliding particles (or scatterers, as the case may be) 
for it to prove the CPR. However, these conditions, under which this construction works to prove the CPR, 
are rather restrictive in the context of NEMD studies of the transport quantities of the systems of physical 
interest, as they exclude a class of very interesting systems where the external force may depend on the 
momenta of the particles — for example, when the system is subjected to an additional external magnetic 
field. Nevertheless, it has allowed us to prove an important theoretical result regarding the CPR, for a gas 
of hard spheres under shear |2^]. It is also possible that this choice of parallel transport can be used for 
numerical computations to achieve higher accuracy for the Lyapunov exponents. 

We end this paper with the observation that while this construction uses explicit dependence on the 
co-ordinate system on a Euclidean manifold, there exists a (co-ordinate independent) differential geometric 
method (for example, in Ref. [16|) that allows one to work easily on non-Euclidean manifolds |23] as well. 
How this present construction can be generalized to such non-Euclidean manifolds is presently unclear, but 
it remains a challenging task for the future. 
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APPENDIX A: 

To simplify the expression of 6p_^_, we break the r.h.s. of Eq. ( pi] ) into several parts. First we notice that 

2 

(I — 2hh) dp_ = y^ Spi- ei+ and (Al) 

i = l 



2 



(I - 2nn) Y,{F* ■ e,_) e,_ 6r - ^(F* • e,+ ) e,+ 5r = - 2 ^ {F* ■ n){e^_ ■ n)ie,^ ■ n) ^^^_ „^^ ^^^^ 



-i — 1 i — 1 i, j — 1 



(eo- • n) 



Next, we observe that the rest of the r.h.s. of Eq. (|3l|), i.e, A • Sr_ describes the effect of the orientation of 
h on 6p_^_. Having denoted the unit vector normal to the surface of the scatterer at C by n , which can be 
related to n by 

h ^ n + dn , (A3) 

{n ■ 6n = 0), it is easily seen from Eqs. (g^26) and ( p8| ) that 



A • Sr_ = - 2 [ (eo- ■ Sn)n + (eo- • n) Sn] . (A4) 

Furthermore, using Sn ~ Sr_/a, where a is the radius of the scatterer. The expression of A • 5r_ can be 
readily simplified as 

a .^-^ (eo-. -n) 

1,3 = 1 
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APPENDIX B: 



Two terms on the r.h.s. of Eq. ( pl[ ) simplify as before, i.e, 



(I - 2N,,,,N,,j,)6P^ = J2 ^P^- 



and 



(Bl) 



ZN-l 



SAf-l 



(I - 2N,^,,N,,j,) Y, (P* ■ e.-) e.- St - J^ i^* ' e»+) ^^+ St 

' {F*-N,,,){h--N,,,){e,--N,,,) 






(eo- • n) 

Following section II, the term A • 5R_ can be expressed as 

A-6R_ = -2[ (eo- • SN^,j,)N,,,, + (eo- • N^,,,)5N,,,,] , 



5Rj- h+ ■ (B2) 



where 



5Ner 



%/2 



0, 0. 



SriiiY , , 



Suiiji , . . . , 



(B3) 



(B4) 



satisfying fiitji ■ Suiiji = 0. This orthogonality condition between hiijr and drii'jr also implies that iV, 



«'j' 



i5iV^ 



0; from which it can be readily shown that 



eo+-[A-dR_] = 0, 



In other words, we have 



3Af-l 



A-5rL= J2{[A■Sli_]■e,+ }e^+. 



(B5) 



(B6) 



The full simplification of the expression [A • SR_] ■ ej+ is carried out below, where it is shown that 

3N-1 



[A-5R_]-e^+ = Y. W,jSRj, 



(B7) 



j = i 



with the property that Wij — Wji (see Eq. ( B18D ). The expression A • SR_, using Eq. (B3), can be written 
in the following form 



A-SR_ 



, , . . . , A,,,. 



A, 



(B8) 



such that all but the i'-th and the j'-th entries on the r.h.s. of Eq.(B8) are zero, while the j'-th entry Aj> is 



related to the i'-th entry Ai' by [using Eq. ( |B4|) ] 



Aj' = - Ai> = [{{pi'- - Pj'-) ■ Sni,j>}ni'j' + {{pi'^ - py-) -hi'j'jSr 



■I- ] 



(B9) 



To obtain an expression for 5niiji analogous to 5n — 5r_/a as used in Appendix A, we need to take a look 
at Figs. 3 and 4. Figure 3 describes, in the laboratory frame, the binary collision process between the i'-th 
and the j'-th sphere on the reference and adjacent trajectories; the thick-lined spheres are on the reference 
trajectory whereas the thin-lined spheres are on the adjacent trajectory. Figure 4 describes the same binary 
collision process in the reference frame of the i'-th sphere (with center C). In Fig. 3, the thick- lined j'-th 
sphere (with center D) on the left depicts the collision situation on the reference trajectory and the thin-lined 
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j'-th sphere (with center E) on the left depicts the coUision situation on the adjacent trajectory. Clearly, in 
Fig. 4, the infinitesimal vector DE is given by 



Sr.^,j, = 5rj,^ - (5ri/_ + [py^ - pe-)ST 
and since the lengths of both the lines CD and CE are 2a (a is the radius of each sphere), we have 

1 



Sn, 



■ Sv ■/ I 
2a '^ 



(BIO) 



(Bll) 



Starting with the expression in Eq. (B6), we proceed to calculate the quantity [A • 5R_] ■ e-ij^. Let us first 
define a 3A^ x 3A^ matrix U composed oi N x N blocks of 3 x 3 matrices, such that, in terms of the block 
indices the only non-zero entries of U are 



V^'^' 



v^r 



-u, 



u, 



I. 



where I is the 3x3 unit matrix. The matrix U has the property that 

U'^ = U and U • N.,,j, == - 2 N,^j, . 
Using Eqs. (B8-B11) and ( B13D , one can then simplify [A • 6R_] ■ ej+ as 

3Af-l 



[A-5R_]-h+ = ^{h^-N,,,) ^ 



2a 



j = i 



(ej_ • U • eo-) - (eo- • U • eo-) 



V2 
2a 



3N-1 r 






(ei+ • U • ej_) - {ei+ ■ U • eo-) 



(e,- 


■N,,,y 


5Ry 


(eo- 


■N,,,,)\ 


. 


^ (ej^-N,,,,) 



{eo--N,'r) 



Equation (B14) can be further simplified by using 

(e,+ • U • eo_) (ej- • iV.-y) = {(e,- • U • eo-) + 4 (e,_ • N^y) (eo- • N,>y)} (e^- • N^y) ., 



(B12) 



(B13) 



SRj- . (B14) 



(B15) 



(ei+ • U • ej-) (eo- • N,,,,) = {(e,- • U • e^--) + 4 (e^- • Ni,j,) (e^- • Ni,j,)} (eo- • N,>j,) 



and Eq. (B13), to obtain 



3N-1 

[A-SR*_]-e^+ = Y. W^jSRj^, 
i = i 



where 



V2a 



(B16) 



(B17) 



(ei- • N,,j,) (cj- • U • eo-) + (e^- • N,>j,) (e,- • U • eo-) - (e,- • U • e^-) (eo- • Ni,j,) 

(eo- -U-eo-) 



(ej--iV,,,')(e,-.iV,,jO 



{eo--N,,,,) 



(B18) 



Equation (B18) is then used in Eq. (B7) 
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p^ 



FIG. 1. Dynamics of the point particle at a collision in two-dimensional projection. 



Adjacent 
trajectory 




Scatterer 



Reference trajectory 



FIG. 2. A two-dimensional projection of the three-dimensional collision dynamics of the reference and adjacent 
trajectories. 
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2a(n././+5n././) 



J- j- 




FIG. 3. Collision between the i'-th and the j'-th sphere on the reference and adjacent trajectories in the labora- 
tory frame. Thick-lined spheres are on the reference trajectory whereas the thin-lined spheres are on the adjacent 
trajectory. 



2a(n,/,/+ 5n././) 



5r._- 5r._+(p;- p,,) 




FIG. 4. Same collisions as in Fig. 3, in the reference frame of the i'-th sphere. 
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